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Abstract 

Background: Chikungunya and dengue viruses emerged in Gabon in 2007, with large outbreaks primarily affecting the 
capital Libreville and several northern towns. Both viruses subsequently spread to the south-east of the country, with new 
outbreaks occurring in 201 0. The mosquito species Aedes albopictus, that was known as a secondary vector for both viruses, 
recently invaded the country and was the primary vector involved in the Gabonese outbreaks. We conducted a 
retrospective study of human sera and mosquitoes collected in Gabon from 2007 to 2010, in order to identify other 
circulating arboviruses. 

Methodology/Principal Findings: Sample collections, including 4312 sera from patients presenting with painful febrile 
disease, and 4665 mosquitoes belonging to 9 species, split into 247 pools (including 137 pools of Aedes albopictus), were 
screened with molecular biology methods. Five human sera and two Aedes albopictus pools, all sampled in an urban setting 
during the 2007 outbreak, were positive for the flavivirus Zika (ZIKV). The ratio of Aedes albopictus pools positive for ZIKV 
was similar to that positive for dengue virus during the concomitant dengue outbreak suggesting similar mosquito 
infection rates and, presumably, underlying a human ZIKV outbreak. ZIKV sequences from the envelope and NS3 genes 
were amplified from a human serum sample. Phylogenetic analysis placed the Gabonese ZIKV at a basal position in the 
African lineage, pointing to ancestral genetic diversification and spread. 

Conclusions/Significance: We provide the first direct evidence of human ZIKV infections in Gabon, and its first occurrence 
in the Asian tiger mosquito, Aedes albopictus. These data reveal an unusual natural life cycle for this virus, occurring in an 
urban environment, and potentially representing a new emerging threat due to this novel association with a highly invasive 
vector whose geographic range is still expanding across the globe. 
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Introduction 

Zika virus (ZIKV) is a mosquito-borne flavivirus phylogeneti- 
cally related to dengue viruses. Following its first isolation in 1947 
from a sentinel monkey placed in the Zika forest in Uganda [1], 
serological surveys and viral isolations (reviewed in [2]) suggested 
that ZIKV (i) ranged widely throughout Africa and Asia, and (ii) 
circulated according to a zoonotic cycle involving non-human 
primates and a broad spectrum of potential mosquito vector 
species. 

In Africa, ZIKV has been isolated from humans in western and 
central countries such as Senegal, Nigeria, Central African 
Republic and Uganda [3-7]. Serological surveys (reviewed in 
[2]) suggested that its geographic range might extend not only to 
other West and Central African countries (Sierra Leone, 
Cameroon, Gabon), but also to eastern (Ethiopia, Kenya, 
Tanzania and Somalia) and northern Africa (Egypt). ZIKV has 
also been isolated from mosquitoes collected in Senegal, Ivory 



Coast, Burkina Faso, Central African Republic and Uganda 
[1,6,8,9]. These mosquitoes mainly belonged to sylvan or rural 
species of the genus Aedes, and more precisely to the Aedimorphus, 
Diceromyia and Stegomyia subgenera. The virus has also been isolated 
in West Africa (Burkina Faso, Senegal and Ivory Coast) [6,9] and 
Asia [10] from Aedes aegypti, a species being considered the main 
ZIKV epidemic vector outside Africa [1 1] . Moreover, Ae. aegypti 
was shown experimentally to be an efficient ZIKV vector [12—14]. 

Despite its apparent broad geographic distribution in Africa and 
Asia, only sporadic cases of human ZIKV infection have been 
reported. This virus received little attention until its sudden 
emergence in Yap Island (Micronesia) in 2007, which involved 
about 5000 persons [15,16], revealing its epidemic capacity. 
Patients develop a mild dengue-like syndrome, including fever, 
headache, rash, arthralgia and conjunctivitis. This clinical similarity 
with other, more commonly diagnosed arboviral infections such as 
chikungunya (CHIKV) and dengue (DENV), might delay the 
diagnosis and/ or lead to underestimation of ZIKV infections. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



1 



February 2014 | Volume 8 | Issue 2 | e2681 



Aedes albopictus Transmission of Zika Virus, Gabon 



Author Summary 

Not previously considered an important human arboviral 
pathogen, the epidemic capacity of Zika virus (ZIKV, a 
dengue-related flavivirus) was revealed by the Micronesia 
outbreak in 2007, which affected about 5000 persons. 
Widely distributed throughout tropical areas of Asia and 
Africa, ZIKV is transmitted by a broad range of mosquito 
species, most of which are sylvatic or rural, Aedes aegypti, 
an anthropophilic and urban species, being considered the 
main ZIKV epidemic vector. In a context of emerging 
arbovirus infections (chikungunya (CHIKV) and dengue 
(DENV)) in Gabon since 2007, we conducted a retrospec- 
tive study to detect other, related viruses. In samples 
collected during the concurrent CHIKV/DENV outbreaks 
that occurred in the capital city in 2007, we detected ZIKV 
in both humans and mosquitoes, and notably the Asian 
mosquito Aedes albopictus that recently invaded the 
country and was the main vector responsible for these 
outbreaks. We found that the Gabonese ZIKV strain 
belonged to the African lineage, and phylogenetic analysis 
suggested ancestral diversification and spread rather than 
recent introduction. These findings, showing for the first 
time epidemic ZIKV activity in an urban environment in 
Central Africa and the presence of ZIKV in the invasive 
mosquito Aedes albopictus, raise the possibility of a new 
emerging threat to human health. 

Here, we report the first direct evidence of ZIKV epidemic 
activity in Central Africa, and its occurrence in an urban 
environment during concomitant CHIKV/DENV outbreaks in 
Libreville, the capital of Gabon, in 2007. We also report the first 
detection of ZIKV in the Asian tiger mosquito, Ae. albopictus. These 
findings, together with the global geographic expansion of this 
invasive species and its increasing importance as epidemic vector 
of arboviruses as exemplified by CHIKV adaptation, suggest that 
the prerequisites for the emergence and global spread of Zika virus 
may soon be satisfied. 

Materials and Methods 

Study 

In 2007 and 2010, Gabon recorded simultaneous outbreaks of 
CHIKV (genus Alphavirus) and DENV (genus Flavivirus) infections. 
The 2007 outbreaks primarily affected Libreville, the capital of 
Gabon, and subsequendy extended northwards to several other 
towns [17], while the 2010 outbreaks occurred in the south-eastern 
provinces [18]. To detect other circulating arboviruses, we 
conducted a retrospective study based on molecular screening of 
4312 sera from symptomatic patients presenting to healthcare 
centers; 24.7% of the samples were obtained during the 2007 
outbreaks, 9.7% during the inter-epidemic period, and 65.5% 
during the 2010 outbreaks (data not shown). We also analyzed a 
collection of 4665 mosquitoes captured during the same period 
and split into 247 pools according to the species, date and 
sampling site (Table 1, see [18] and [19] for the details of the 
methodology used for mosquito trapping). 

Ethics statement 

The Centre International de Recherches Medicales de France- 
ville (CIRMF) and the Gabonese Ministry of Health cooperated in 
the 2007 and 2010 outbreak response and management, that 
included blood sampling for laboratory diagnostic and epidemi- 
ological survey. The study was approved by our Institutional 
review board (Conseil scientifique du CIRMF). 



Symptomatic patients presented to health care centers for 
medical examination. All patients were informed that blood 
sampling was required for laboratory diagnosis of suspected acute 
infections, such as malaria, dengue or chikungunya fever. During 
the two outbreaks, given the urgency of diagnosis, only oral 
consent was obtained for blood sampling and was approved by the 
institutional review board. However during the active surveillance 
study that was performed between the two outbreaks (described in 
reference [18]), written consent could be obtained. 

Virus identification and characterization 

Primary molecular screening was based on hemi-nested reverse- 
transcription PCR (hnRT-PCR) with the generic primers PF1S/ 
PF2Rbis/PF3S targeting highly conserved motifs in the flavivirus 
polymerase (NS5) gene (280-bp) [20]. Yellow fever virus RNA 
(vaccinal strain 17D) was used as a positive control. A second 
screening was performed with a ZIKV-specific real-time PCR 
method using the primers-probe system ZIKV-1086/ZIKV- 
1162c/ZIKV-1107-FAM [16], also targeting a short sequence 
(160 bp) of the NS5 gene. 

Virus isolation was attempted on the Vero and C6/36 cell lines 
but was unsuccessful, presumably because of low viral titers 
(despite two patients presenting only 1 and 4 days after symptom 
onset), and unsuitable initial storage conditions. To further 
characterize the Gabonese ZIKV strains, partial envelope (E) 
(841 bp) and NS3 (772 bp) gene sequences were amplified by 
conventional nested RT-PCR with specific primers derived from 
published ZIKV sequences. The primer pairs targeting the E gene 
were ZIK-ES1 (TGGGGAAAYGGDTGTGGACTYTTTGG)/ 
ZIK-ER1 (CCYCCRACTGATCCRAARTCCCA) and ZIK- 
ES2 (GGGAGYYTGGTGACATGYGCYAAGTT)/ZIK-ER2 
(CCRATGGTGCTRCCACTCCTRTGCCA). The primer pairs 
for NS3 amplification were ZIK-NS3FS (GGRGTCTTCCACA- 
CYATGTGGCACGTYACA)/ZIK-NS3FR (TTCCTGCCTAT- 
RCGYCCYCTCCTCTGRGCAGC) and ZIK-X1 (AGAGTGA- 
TAGGACTCTATGG)/ZIK-X2 (GTTGGCRCCCATCTCT- 
GARATGTCAGT). 

Phylogenetic analysis 

The E and NS3 sequences obtained from one Gabonese patient 
were concatenated and analyzed using a set of previously 
published ZIKV sequences. Phylogenetic relationships were 
reconstructed with the maximum likelihood algorithm implement- 
ed in PhyML [21] (available at http://www.atgc-montpellier.fr/ 
phyml/) with best of NNI (Nearest Neighbor Interchange) and 
SPR (Subtree Pruning and Regrafting) criteria for tree topology 
searching, and the GTR model of nucleotide substitutions. The 
Gamma distribution of rate heterogeneity was set to 4 categories, 
with a proportion of invariable sites and an alpha parameter 
estimated from the dataset. Branch support was assessed from 1 00 
bootstrap replicates. Tree reconstructions were also performed by 
Bayesian inference with MrBayes v3.2 [22] under the GTR+I+G 
model of nucleotide substitutions, and with the distance neighbor- 
joining method [23] implemented in MEGA5 [24] with 
confidence levels estimated for 1000 replicates. To test for 
phylogenetic discrepancies, tree reconstructions were also per- 
formed independently from the envelope dataset and the NS3 
dataset with PhyML according to the parameters described above. 

The resulting trees were visualized with the FigTree software 
(Available at: http:/ /tree. bio. ed.ac.uk/software/figtree/), and 
rooted on midpoint for clarity. The Genbank accession numbers 
for the Gabonese ZIKV strain are KF270886 (envelope) and 
KF270887 (NS3). 
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Table 1. Mosquito collections screened for Zika virus. 



Libreville 2007 












Franceville 
2010 


Total 




Species Pools 


Mos. 


Id. (No.) 


ZIKV 


CHIKV 


DENV 


Libreville 
suburb 


Pools Mos. 


Pools (%) Mos. 


Aedes albopictus 91 


2130 


T64 (21) 


+ 


+ 




Nzeng-Ayong 


46 571 


137 (55.4) 


2701 






T713 (25) 


+ 


+ 




Alenkiri 












T707 (25) 






+ 


Alenkiri 












T717 (25) 






+ 


Alenkiri 












T723 (25) 






+ 


Alenkiri 












T724 (6) 




+ 




Alenkiri 












T21 (25) 




+ 




Avorembam 












T22 (25) 




+ 




Avorembam 












T280 (1) 




+ 




Bel-Air 








Aedes aegypti* 40 


853 












5 28 


45 (18.2) 


881 


Aedes simpsoni 10 
complex 


52 












5 36 


15 (6.1) 


88 


Anopheles gambiae* 8 


72 














8 (3.2) 


72 


Mansonia africana 6 


86 














6 (2.4) 


86 


Mansonia uniformis* 4 


99 














4 (1.6) 


99 


Culex quinquefasciatus 29 


690 














29 (11.7) 


690 


Culex spp. 1 


22 














1 (0.4) 


22 


Eretmapodites 
quinquevittatus* 














2 26 


2 (0.8) 


26 


Total 189 


4004 












58 661 


247 


4665 



*Species in which Zika virus has previously been detected. 

(%) The percentage of each mosquito species in the collection is indicated in brackets. Mos.: Number of mosquitoes included in a pool. Id. (No.): Mosquito pool positive 
for ZIKV, CHIKV or DENV, followed by the total number of included mosquitoes in the pool indicated in brackets. 
doi:1 0.1 371 /journal.pntd.0002681 .t001 



Results 

Molecular screening 

The NS5 PCR products were sequenced, resulting in the first 
ZIKV RNA detection in a human sample (Cocobeach) and in two 
Ae. albopictus pools (Libreville) collected during the 2007 outbreaks. 
Real-time PCR was then performed, leading to the detection of 
four additional positive human samples, collected in 2007 in four 
suburbs of Libreville (Diba-Diba, Nzeng-Ayong, PK8, PK9) 
(Figure 1). No ZIKV was detected during the inter-epidemic 
period or during the 2010 outbreaks. 

Clinical description 

Clinical information was available for only one ZIKV-positive 
patient, who had mild arthralgia, subjective fever, headache, rash, 
mild asthenia, myalgia, diarrhea and vomiting. No information 
was available on this patient's outcome. Cycle threshold values for 
human blood samples were high (>37 cycles), suggesting low viral 
loads (data not shown). 

Vector involvement 

Aedes albopictus was the predominant species collected, account- 
ing for 55.4% of the mosquito pools, while Aedes aegypti accounted 
for 18.2% (Table 1). The other mosquito species consisted of 
members of the Aedes simpsoni complex, Anopheles gambiae, Mansonia 
africana, Mansonia uniformis, Culex quinquefasciatus, Eretmapodites 
quinquevittatus and unidentified Culex species. Positive mosquito 
pools were captured from two suburbs (Nzeng-Ayong and 



Alenkiri) where Aedes albopictus was the predominant species 
(Figure 1, Table 1). 

Sequences analysis 

As isolation on the Vero and C6/36 cell lines failed, the 
Gabonese ZIKV strain was further characterized by partial 
sequencing of the E and NS3 genes. Phylogenetic analysis was 
performed on concatenated E and NS3 sequences from one 
Cocobeach serum sample. The resulting tree topology (Figure 2) 
was similar to that previously obtained from the complete coding 
sequences, corroborating Asian and African distinct lineages [2]. 
The African lineage was further split into two groups, one 
containing the genetic variants of the MR766 strain (Uganda, 
1947) and the second including West African strains (Nigeria, 
1968; Senegal, 1984) and the new ZIKV sequence from Gabon, at 
a basal position. Phylogenetic trees derived from the E and NS3 
partial sequences resulted in a similar topology, apart from the 
weakly supported branching pattern for the MR766 variant 
DQ859059, oscillating between the two African sister groups 
(Supporting Figure SI). The deletions in potential glycosylation 
sites previously reported for the Nigerian ZIKV strain and two 
variants of the Ugandan strain MR766 (sequences AY632535 and 
DQ859059) [2] were absent from the Gabonese ZIKV sequence. 

Discussion 

Evidence of human ZIKV infections in Central Africa is limited to 
one isolate from RCA in 1991 [6] and two serological surveys 
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Figure 1. Geographic distribution of Zika and chikungunya and/or dengue viruses infections in Gabon in 2007. The left-hand panel 
indicates Gabonese CHIKV and/or DENV cases in green circles and ZIKV cases in purple circles. The right-hand panel shows the location of Libreville 
suburbs where ZIKV-positive human sera (H) and mosquito pools (M) were detected. 
doi:10.1371/journal.pntd.0002681.g001 



performed 50 years ago in Gabon [25,26]. No report of human cycle. We provide here the first direct evidence of human ZIKV 
ZIKV infections was made in other countries of the Congo basin infections in Gabon, as well as its occurrence in an urban transmission 
forest block, despite probable circulation through a sylvan natural cycle, and the probable role of Ae. albopictus as an epidemic vector. 



100/100/100 



Malaisia 1966 (HQ234499) 

r— Micronesia 2008 (EU545988) 
I Cambodia 2010 (JN860885) 



100/100/100 



66/87/65 



100/100/100 



100/100/100 



substitution/site 
0,03 



100/100/100 
92/91/97 



Uganda 1947 (DQ859059) 
Uganda 1947 (AY632535) 
Uganda 1947 (AY234498) 
— Gabon 2007 



I Nigeria 1968 (HQ234500) 



-Senegal 1984 (HQ234501) 



Figure 2. Phylogenetic relationships between concatenated sequences of the Zika virus envelope and NS3 genes. The tree was 
constructed with the maximum likelihood algorithm implemented in PhyML and rooted on midpoint. Bootstrap values are shown at the respective 
nodes, followed by bootstrap values resulting from NJ analysis and, finally, the posterior probability resulting from Bayesian analysis. The scale bar 
indicates the number of substitutions per site. The GenBank accession numbers for the 2007 Gabonese ZIKV isolate are KF270886 (envelope) and 
KF270887 (NS3). 

doi:1 0.1 371 /journal.pntd.0002681 ,g002 
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Our phylogenetic results are in agreement with the tree topology 
previously obtained with complete coding sequences of ZIKV strains, 
showing an African lineage and an Asian lineage [2] . The branching 
pattern obtained here suggests that ZIKV emergence in Gabon did 
not result from strain importation but rather from the diversification 
and spread of an ancestral strain belonging to the African lineage. 
The identification of ZIKV in two different localities of Gabon 
(Cocobeach and Libreville) suggests that the virus was widespread 
rather than restricted to a single epidemic focus. The simultaneous 
occurrence of human and mosquito infections in Libreville also 
suggests that the virus circulated in 2007 in an epidemic cycle rather 
than as isolated cases introduced from sylvan cycles. 

Of note, ZIKV transmission occurred here in a previously 
undocumented urban cycle, supporting the potential for urbani- 
zation suggested in 2010 by Weaver and Reisen [27]. While some 
mosquito species (including Ae. aegypti) previously found to be 
associated with ZIKV, were captured and tested here, only Ae. 
albopictus pools were positive for this virus. Moreover, this species 
largely outnumbered Ae. aegypti in the suburbs of Libreville where 
human cases were detected, suggesting that Ae. albopictus played a 
major role in ZIKV transmission in Libreville. 

The ratio of ZIKV-positive Ae. albopictus pools is similar to that 
reported for DENV-positive pools, suggesting that these two 
viruses infect similar proportions of mosquitoes. The small number 
of recorded human ZIKV cases, by comparison with DENV cases, 
may be due to the occurrence of subclinical forms of ZIKV 
infections that did not required medical attention. Thus, an 
underlying ZIKV epidemic transmission might have been masked 
by concomitant CHIKV/DENV outbreaks. 

The natural histories of CHIKV and ZIKV display several 
similarities. Before the large Indian Ocean outbreaks in 2005-2007, 
chikungunya fever was a neglected arboviral disease. Both viruses 
are phylogenetically closely related to African viruses [28-30] 
suggesting they probably originated in Africa, where they circulated 
in an enzootic sylvan cycle involving non-human primates and a 
wide variety of mosquito species, human outbreaks presumably 
being mediated by Ae. aegypti [5,31]. In Asia, both viruses are 
thought to circulate mainly in a human-mosquito cycle involving^. 
aegypti [11,14,31]. Together with the recent Yap Island outbreak, 
this prompted some researchers to re-examine the susceptibility of 
Ae. aegypti to ZIKV infection [14]. However, it must be noted that 
the vector of the Yap Island outbreak was not definitely identified 
since the predominant potential vector species Aedes hensilli remained 
negative [15], and that ZIKV has been isolated only once from Ae. 
aegypti in Asia [ 1 0] , so that its vector status in natura is not confirmed. 
Additionally, a ZIKV enzootic transmission cycle involving non- 
human primates in Asia and sylvatic vectors cannot be ruled out as 
suggested by serologic studies carried on orangutang [32,33]. 
Finally both CHIKV and ZIKV have shown their ability to adapt to 
a new vector, Ae. albopictus, upon introduction in an environment 
where their primary vector was outnumbered. This mosquito 
species being native to South-East Asia, our findings may help to 
explain human ZIKV transmission in Asia. 

Aedes albopictus was first introduced in Africa in 1991 [34] and 
detected in Gabon in 2007, where its invasion likely contributed to 
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